Method for fault location on series compensated power transmission lines with two-end unsynchronized measurement

ABSTRACT

A method for fault location in series compensated power transmission lines using two-end unsynchronized voltage and current measurements in stations A and B. The method includes performing a subroutine I for determining whether the fault occurred in a line section LA between the station A and a SC&amp;MOV and the distance d A  to the fault from station A. Subroutine II determines whether the fault occurred in a line section LB between the station B and the SC&amp;MOV and the distance d B  to the fault from station B. The synchronization angle e jδ     A    is determined in subroutine I, and the synchronization angle e jδ     B    is determined in subroutine II. Then, the distance to fault d A  or d B  is determined. Finally, an equivalent impedance of compensating banks at fault stage  Z   SC1     —     ph   SUB     —     A  and  Z   SC1     —     ph   SUB     —     B , and pre-fault stage  Z   SC1     —     pre  is calculated in order to determine whether distance d A  or d B  is the final result.

BACKGROUND OF THE INVENTION

The present invention is concerned with a method for fault location on series compensated power transmission lines with two-end unsynchronized measurement, finding an application in the electrical power sector and for overhead transmission lines uncompensated and with series compensation.

The location fault technique on series compensated power transmission lines is known from EP1051632B1. In this solution it has been proposed to calculate distance to fault with utilizing local measurements from one terminal of the series compensated transmission line where fault locator is placed. This method makes use of transforming the parallel connection of series capacitor and metal oxide varistor (MOV) into a model—equivalent impedance in the form of a series connection of resistance and reactance, both fundamental current magnitude dependent. This equivalent impedance is used to calculate the distance to fault. This method requires the knowledge of the parameters of series capacitor banks and MOV elements as well.

A method of reach measurement for distance relay and fault locators for series—compensated transmission lines is know from U.S. Pat. No. 6,336,059 B1. In the method according to that invention the fault location estimation is done with using local information only. In the method a number of line currents samples are measured, where such samples are representative of values of line current waveform at successive instants of time on the series compensation line. Capacitor voltage values are computed based on measured currents in accordance with an equation which takes into account the non-linear protection device—metal oxide varistor parallel to the installed series capacitor. The distance to fault is calculated with taking into account the calculated earlier capacitor voltage values.

In US patent 2006/0142964 A1 it is presented the fault location method with using two end currents and voltages measurement. The invention is characterised in that after determination of current transformer saturation at first end of the line A or second end B, the distance to fault is calculated with using voltage from both ends of the line and currents only from this side of the line where current transformers are not saturated. Calculation of distance to fault is done based on generalized fault loop model. In this model the total fault current is determined with use of lumped line model and, one end current signals and two ends voltage signals. However the method according to the cited invention is not designed for locating faults on series compensated lines.

Another method of fault location on series compensation lines with using two-end measurements is known from the article Chi-Shan Yu, Chih-Wen Liu, Sun-Li Yu, and Joe-Air Jiang, “A New PMU-Based Fault Location Algorithm for Series Compensated Lines”, IEEE Transactions on Power Delivery, VOL. 17, NO. 1, JANUARY 2002, pp. 33-46. This method uses voltage and current signals measured synchronously at both ends of a single series compensated line. Thus, differs by the kind of measurement (synchronous measurement) and type of a line (single series-compensated line) from the invented approach, which is designed for use of more general case of asynchronous measurements and in application to both single and double-circuit uncompensated and series-compensated lines. Moreover, the cited approach does not provide an analytic formula for the sought distance to fault but is based on scanning along the whole line section for finding the fault position at which the determined fault voltage and total fault current are in phase (this is so since the fault path has the resistive character). For each checked point on the line, the fault voltage and total fault current have to be determined from the symmetrical components of current and voltage. This imposes high computational burden required for performing the fault location. In contrast, the invented approach delivers the compact formula for the sought distance to fault, derived with strict considering the distributed parameter model of the line.

From the article Claude Fecteau “Accurate Fault Location Algorithm for Series Compensated Lines Using Two-Terminal Unsynchronized Measurements and Hydro-Quebec's Field Experience”, presented to the 33-rd Annual Western Protective Relay Conference, Spokane, Wash., Oct. 17-10, 2006 Conference there is known the solution that is to a certain extent related with the invented approach. One can point out that the developed approach is more general (single and double-circuit lines, uncompensated and series-compensated lines). The cited approach does not provide an analytic formula for the sought distance to fault. Instead, the distance to fault is determined iteratively by minimizing the objective function for the reactance of the fault impedance. This imposes higher computational burden than the invented approach.

Moreover, the invented approach offers additionally determination of the synchronization angle, in case of single phase-to-ground faults and phase-to-phase faults (majority of the faults in practice), with use of the post-fault quantities. This is advantageous from the assuring precise synchronization point of view.

SUMMARY OF THE INVENTION

The essence of the inventive method for fault location in series compensated power transmission lines with two-end unsynchronized measurement wherein a distance to fault and fault resistance are determined by means of voltages and currents values measurement in the stations A and B, before and after occurrence of the fault, consist of the following:

-   -   if the fault occurred in the line section LA between station A         and series capacitor and metal oxide varistor device SC&MOV, the         distance for fault location d_(A) is determined by         -   determination of a synchronization angle δA in term e^(jδ)             ^(A) for known type of the fault which is phase-to-ground             faults or phase-to-phase faults, or which is not the             phase-to-ground faults or not the phase-to-phase faults,         -   then taking into account the distributed parameter line             model and solving the formula:             d _(A) =p _(SC) d _(LA),         -   where:         -   p_(SC)—denotes the distance (p.u.) per unit of total line             length l at which the compensation bank SC&MOV is installed             far away from the station A,         -   d_(LA)—denotes an hypothetical distance to fault, expressed             in relation to length p_(SC)l and which value is received             using known iterative method by solving the fault loop             formula:

${{\sum\limits_{i = 0}^{2}{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Ai}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Ai}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}}} \right)} \cdot {\mathbb{e}}^{j\;\delta_{A}}}} - {R_{FA}{\sum\limits_{i = 1}^{2}\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LA}} \right)}p_{SC}l} \right)}}}} = 0$ for single line,

${{\sum\limits_{i = 0}^{2}{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Ai}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Ai}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}}} \right)} \cdot {\mathbb{e}}^{j\;\delta_{A}}}} - {R_{FA}{\sum\limits_{i = 1}^{2}\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LA}} \right)}p_{SC}l} \right)}}} - {{\underset{\_}{a}}_{m\; 0}d_{LA}p_{SC}{\underset{\_}{Z}}_{0\; m}{{\underset{\_}{I}}_{{{Aparal}\_}0} \cdot {\mathbb{e}}^{j\;\delta_{A}}}}} = 0$ for double-circuit line, in which: R_(FA)—denotes the unknown fault resistance, a ₁, a ₂, a ₀, a _(m0)—denotes the weighting coefficients dependent on fault type, M _(i)—denote numerical components for positive and negative sequences, a _(F1), a _(F2)—denotes the share coefficients dependent on fault type, Z _(0m)—mutual coupling impedance for zero sequence, I _(Aparal) _(—) ₀—denotes zero sequence current from the station A measured in healthy parallel line,

-   -   then an equivalent impedance of compensating bank at fault stage         (Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A)) is calculated with using         the post fault of measured voltage and currents and an         equivalent impedance of compensation bank at pre-fault (Z_(SC1)         _(—) _(pre)) is calculated in order to determine the final         result d_(A) which is chosen in the following steps:     -   first, it is check whether the fault resistances R_(FA) is of         positive value and if not, then such subroutine is rejected,     -   next, the features “real” and:“imag” of the estimated equivalent         circuit of compensation banks SCs&MOVS are checked and if each         of it is satisfied the following dependence:         real(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A))>0,         imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A))<0,         and imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A))>imag(Z _(SC1)         _(—) _(pre))

then

-   -   that the fault has occurred between the station A and the SC&MOV         bank, and the result d_(A) is a valid result.

The essence of the inventive method for fault location in series compensated power transmission lines with two-end unsynchronized measurement wherein a distance to fault and fault resistance are determined by means of voltages and currents values measurement in the stations A and B, before and after occurrence of the fault, consist of the following:

-   -   if the fault occurred in the line section LB between station B         and series capacitor and metal oxide varistor device SC&MOV, the         distance d_(B) for fault location id determined by         -   determination of a synchronization angle δ_(B) in term             e^(jδ) ^(B) for known type of the fault which is             phase-to-ground faults or phase-to-phase faults, or which is             not the phase-to-ground faults or not the phase-to-phase             faults,         -   then taking into account the distributed parameter line             model and solving the formula:             d _(B)=(1−P _(SC))d _(LB)         -   where:         -   (1−p_(SC))—denotes the distance (p.u.) per unit of total             line length l, at which the compensation bank SC&MOV is             installed far away from the station B,             -   d_(LB)—denotes an hypothetical distance to fault,                 expressed in relation to length (1−p_(SC)l) and which                 value is received using known iterative method by                 solving the fault loop formula:

${{\sum\limits_{i = 0}^{2}{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Bi}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Bi}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}}} \right)} \cdot {\mathbb{e}}^{j\;\delta_{B}}}} - {R_{FB}{\sum\limits_{i = 1}^{2}\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LB}} \right)}\left( {1 - p_{SC}} \right)l} \right)}}}} = 0$ for single line,

${{\sum\limits_{i = 0}^{2}{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Bi}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Bi}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}}} \right)} \cdot {\mathbb{e}}^{j\;\delta_{B}}}} - {R_{FB}{\sum\limits_{i = 1}^{2}\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LB}} \right)}\left( {1 - p_{SC}} \right)l} \right)}}} - {{\underset{\_}{a}}_{m\; 0}{d_{LB}\left( {1 - p_{SC}} \right)}{\underset{\_}{Z}}_{0\; m}{{\underset{\_}{I}}_{{Bparal}\;\_ 0} \cdot {\mathbb{e}}^{j\;\delta_{B}}}}} = 0$ for double-circuit line where: R_(FB)—denotes the unknown fault resistance, M _(i)—denote numerical coefficients for positive and negative sequences, a ₁, a ₂, a ₀, a _(m0)—denotes the weighting coefficients dependent on fault type, a _(F1), a _(F2)—denotes the share coefficients dependent on fault type. Z _(0m)—mutual coupling impedance for zero sequence, I _(Bparal) _(—) ₀—denotes zero sequence current from the station B measured in healthy parallel line,

-   -   then an equivalent impedance of compensating bank at fault stage         (Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B)) is calculated with using         the post fault of measured voltage and currents and an         equivalent impedance of compensation bank at pre-fault (Z_(SC1)         _(—) _(pre)) is calculated in order to determine the final         result d_(B) which is chosen in the following steps:     -   first, it is check whether the fault resistances R_(FB) are of         positive value and if not, then such subroutine is rejected,     -   next, the features “real” and:“imag” of the estimated equivalent         circuit of compensation banks SCs&MOVS are checked and if it is         satisfied the following dependence:

real(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B))>0

imag(Z_(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B))<0

and imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B))>imag(Z _(SC1) _(—) _(pre))

then that the fault has occurred between the station B and the SC&MOV bank, and the result d_(B) is a valid result.

Preferably the synchronization angle δ_(A) in term e^(jδ) ^(A) is determined for phase-to-ground faults or phase-to-phase faults from the formula:

$\left\lbrack {\mathbb{e}}^{j\;\delta_{A}} \right\rbrack_{{{ph} - g},{{ph} - {ph}}} = \frac{{{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{I}}_{{SC}\; 2\_\; B}} - {{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{I}}_{{SC}\; 1\_\; B}}}{{{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{N}}_{A\; 1}} - {{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{N}}_{A\; 2}}}$

-   -   where:     -   subscript ph-g denotes phase to ground faults: a-g, b-g, c-g,         and subscript ph-ph denotes phase to phase faults a-b, b-c, c-a,     -   a _(F2) ^(I-SET), a _(F1) ^(II-SET)—denote the coefficients         dependent on fault type,         I _(SC1) _(—) _(B)—denotes the positive sequence current from         the station B after analytical transfer towards series capacitor         banks SC&MOV,         I _(SC2) _(—) _(B)—denotes the negative sequence current from         the station B after analytical transfer towards series capacitor         banks SC&MOV,         N _(A1), N _(A2)—denote the coefficients dependent on positive         and negative sequence currents and voltage taken at station A         respectively as well as parameters of the line.         Preferably synchronization angle δ_(A) in term e^(jδ) ^(A) is         determined for the faults nor being the phase-to-ground faults         or phase-to-phase faults from the formula:

${\mathbb{e}}^{j\;\delta_{A}} = \frac{{\left( {1 + {{\underset{\_}{Z}}_{1\; B}{\underset{\_}{Y}}_{1\; B}}} \right)\left( {{{\underset{\_}{Y}}_{1\; B}{\underset{\_}{V}}_{B\; 1}^{pre}} - I_{B\; 1}^{pre}} \right)} + {{\underset{\_}{Y}}_{1\; B}{\underset{\_}{V}}_{B\; 1}^{pre}}}{{\left( {1 + {{\underset{\_}{Z}}_{1\; A}{\underset{\_}{Y}}_{1\; A}}} \right)\left( {{\underset{\_}{I}}_{A\; 1}^{pre} - {{\underset{\_}{Y}}_{A\; 1}^{pre}V_{A\; 1}^{pre}}} \right)} - {{\underset{\_}{Y}}_{1\; A}{\underset{\_}{V}}_{A\; 1}^{pre}}}$ where: V _(A1) ^(pre), V _(B1) ^(pre)—denote pre-fault positive sequence quantities (superscript ‘pre’ and subscript ‘1’) of voltage measured at the ends A and B, respectively, I _(A1) ^(pre), I _(B1) ^(pre)—denote pre-fault positive sequence quantities (superscript ‘pre’ and subscript ‘1’) of current measured of the faulted line at the ends A and B, respectively, Z _(1A), Y _(1A)—denote impedance and admittance of line section LA, Z _(1B), Y _(1B)—denote impedance and admittance of line section LB, LA—is the part of the line between line station A and series capacitor banks LB—is the part of the line between line station B and series capacitor banks. Preferably synchronization angle δ_(B) in term e^(jδ) ^(B) is determined for phase-to-ground faults or phase-to-phase faults from the formula:

$\left\lbrack {\mathbb{e}}^{{j\delta}_{B}} \right\rbrack_{{{ph} - g},{{ph} - {ph}}} = \frac{{{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{I}}_{{SC}\; 2{\_ A}}} - {{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{I}}_{{SC}\; 1{\_ A}}}}{{{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{N}}_{B\; 1}} - {{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{N}}_{B\; 2}}}$

where:

a _(F2) ^(I-SET), a _(F1) ^(II-SET)—denote the coefficients dependent of fault type,

I _(SC1) _(—) _(A)—denotes positive sequence current from the station A after analytical transfer towards series capacitor banks SC&MOV,

I _(SC2) _(—) _(A)—denotes negative sequence current from the station A after analytical transfer towards series capacitor banks SC&MOV,

N _(B1), N _(B2)—denote the coefficients dependent on positive and negative sequence currents and voltage taken at station B respectively as well as parameters of the line.

Preferably the synchronization angle δ_(B) in term e^(jδ) ^(B) is determined for the faults nor being the phase-to-ground faults or phase-to-phase faults from the formula:

${\mathbb{e}}^{- {j\delta}_{B}} = \frac{{\left( {1 + {{\underset{\_}{Z}}_{1B}{\underset{\_}{Y}}_{1B}}} \right)\left( {{{\underset{\_}{Y}}_{1B}{\underset{\_}{V}}_{B\; 1}^{pre}} - {\underset{\_}{I}}_{B\; 1}^{pre}} \right)} + {{\underset{\_}{Y}}_{1B}{\underset{\_}{V}}_{B\; 1}^{pre}}}{{\left( {1 + {{\underset{\_}{Z}}_{1A}{\underset{\_}{Y}}_{1A}}} \right)\left( {{\underset{\_}{I}}_{A\; 1}^{pre} - {{\underset{\_}{Y}}_{1A}{\underset{\_}{V}}_{A\; 1}^{pre}}} \right)} - {{\underset{\_}{Y}}_{1A}{\underset{\_}{V}}_{A\; 1}^{pre}}}$ where: V _(A1) ^(pre), V _(B1) ^(pre)—denote pre-fault positive sequence quantities (superscript ‘pre’ and subscript ‘1’) of voltage measured at the ends A and B, respectively, I _(A1) ^(pre), I _(B1) ^(pre)—denote pre-fault positive sequence quantities (superscript ‘pre’ and subscript ‘1’) of current measured of the faulted line at the ends A and B, respectively, Z _(1A), Y _(1A)—denote impedance and admittance of line section LA, Z _(1B), Y _(1B)—denote impedance and admittance of line section LB. LA—is the part of the line between line station A and series capacitor banks LB—is the part of the between line station B and series capacitor banks.

The advantage of the method according to the invention is the overcoming of all the limitations and drawbacks of the known methods that means that parameters and status of the series capacitor bank are not required to be known since they are not utilized for determination of the distance to fault and fault resistance, i.e. only the position at which the compensating banks are installed has to be known. Additionally the reactance of series capacitors is measured using the pre-fault measurements. The synchronization angle is determined using the pre-fault measurements or alternatively, for single phase-to-ground faults and phase-to-phase faults which are the most common faults with using the post-fault measurements, and with strict considering the distributed parameter line model, what allows accurate analytical synchronization of the measurements acquired asynchronously. Then, the distances to fault are calculated by introducing the compensation for shunt capacitances of the line with use of the distributed parameter line model. Selection of the valid subroutine is done based on multi-criteria procedure, what allows for reliable valid results identification within very wide range of fault conditions, even for very high fault resistances.

BRIEF DESCRIPTION OF THE DRAWINGS

The method according to the present invention is explained in an embodiment shown in the drawing, where

FIG. 1 presents a general diagram of the electric network for the implementation of the inventive method,

FIG. 2—an equivalent circuit diagram of series-compensated line for pre-fault positive sequence considered for determining the synchronization angle,

FIG. 3—for fault FA an equivalent circuit diagram of series-compensated line for positive sequence—the part containing the SC&MOV and the section between the SC&MOV and the station B,

FIG. 4—for fault FA an equivalent circuit diagram of series-compensated line for positive sequence—the part containing the SC&MOV and the section between the SC&MOV and the station A,

FIG. 5—for fault FA an equivalent circuit diagram of series-compensated line for negative sequence—the part containing the SC&MOV and the section between the SC&MOV and the station B,

FIG. 6—for fault FA an equivalent circuit diagram of series-compensated line for negative sequence—the part containing the SC&MOV and the section between the SC&MOV and the station A,

FIG. 7—for fault FB an equivalent circuit diagram of series-compensated line for positive sequence—the part containing the SC&MOV and the section between the SC&MOV and the station A,

FIG. 8—for fault FB an equivalent circuit diagram of series-compensated line for positive sequence—the part containing the SC&MOV and the section between the SC&MOV and the station B,

FIG. 9—for fault FB an equivalent circuit diagram of series-compensated line for negative sequence—the part containing the SC&MOV and the section between the SC&MOV and the station A,

FIG. 10—for fault FB an equivalent circuit diagram of series-compensated line for negative sequence—the part containing the SC&MOV and the section between the SC&MOV and the station B,

FIG. 11—an equivalent circuit diagram of series—compensated line for the pre-fault positive sequence, considering for measuring the compensating bank impedance,

FIG. 12—a flow chart of the example of a fault location algorithm according to the present invention.

DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS

The transmission system shown in FIG. 1 consists of series compensated transmission line connecting two electric power stations A, and B. The installed series capacitor (SC-series capacitor) includes a metal oxide varistor (MOV—metal oxide varistor) connected in parallel (SC&MOV) which limits the voltage across the capacitor in a pre-defined manner. Station A is located at the beginning of the line, station B at the end of this line. In station A there is a fault locator FL. The fault locator can be also located at station B or at the separate place, what is not shown at the FIG. 1. Fault location is done using distributed parameters line model, models of faults and fault loops for symmetrical components positive, negative and zero, and different types of faults, by applying suitable share coefficients determining the relation between the symmetrical components of the total fault current when voltage drop across the fault resistance is estimated, defined as a _(F1), a _(F2), a _(F0) and weight coefficients a ₁, a ₂, a ₀, a _(m0) defining the share of individual components in the total model of the fault loop.

Unsynchronised measurements of three-phase currents from stations A, B, and of three-phase voltages from station A and B are supplied to the fault locator FL. Fault locator is supplied with information about the type of the fault and the time of its occurrence. Then two hypothetical fault points are assumed and the distances d_(A), and d_(B) between the beginning of the given line section A and B and the fault points (FA, FB) are calculated in order to determinate valid results, i.e. consistent with the real fault case, on the following assumptions:

-   -   calculation of synchronization angle δ_(A), the distance d_(A)         from the beginning of the line in the station A to the fault         point FA, fault resistance R_(FA) and impedance of compensation         bank at faults state Z _(SC1) _(—) _(ph) ^(SUB) ^(—A) ,         impedance of compensation bank at pre-fault state Z _(SC1) _(—)         _(pre), assuming that the fault occurred in line section         LA—steps 103 a, 104 a, 105 a, 106 a, where LA is the part of the         line where the fault occurred between line station A and         capacitor bank,     -   calculation of synchronization angle δ_(B), the distance d_(B)         from the beginning of the line in the station B to the fault         point FB, fault resistance R_(FB) and impedance of compensation         bank at faults stage Z _(SC1) _(—) _(ph) ^(SUB) ^(—A) ,         impedance of compensation bank at pre-fault state Z _(SC1) _(—)         _(pre), assuming that the fault occurred in line section LB,         where LB is the part of the line where the fault occurred         between line station B and capacitor bank,—steps 103 b, 104 b,         105 b, 106 b.         Then among those two distances to fault d_(A) and d_(B) the         valid results is selected—step 107.         Step 101. In stations A and B, current and voltage input signals         from individual phases for fault conditions I_(A) _(—) _(a),         I_(A) _(—) _(b), I_(A) _(—) _(c), V_(A) _(—) _(a), V_(A) _(—)         _(b), V_(A) _(—) _(c), I_(B) _(—) _(a), I_(B) _(—) _(b), I_(B)         _(—) _(c), V_(B) _(—) _(a), V_(B) _(—) _(b), V_(B) _(—) _(c) and         for pre-fault conditions I_(A) _(—) _(pre) _(—) _(a), I_(A) _(—)         _(pre) _(—) _(b), I_(A) _(—) _(pre) _(—) _(c), V_(A) _(—) _(pre)         _(—) _(a), V_(A) _(—) _(pre) _(—) _(b), V_(A) _(—) _(pre) _(—)         _(c), I_(B) _(—) _(pre) _(—) _(a), I_(B) _(—) _(pre) _(—) _(b),         I_(B) _(—) _(pre) _(—) _(c), V_(B) _(—) _(pre) _(—) _(a), V_(B)         _(—) _(pre) _(—) _(b), V_(B) _(—) _(pre) _(—) _(c) and in case         of double-circuit line currents from the healthy parallel line         I_(Aparal) _(—) _(a), I_(Aparal) _(—b) , I_(Aparal) _(—) _(c),         I_(Bparal) _(—) _(a), I_(Bparal) _(—) _(b), I_(Bparal) _(—) _(c)         or only zero sequence currents I_(Aparal) _(—) ₀, I_(Bparal)         _(—) ₀ are measured.         Step 102. In this step, the symmetrical components of the phase         currents measured in stations A and B, and of phase voltages         measured in stations A and B for pre-fault conditions and fault         conditions are calculated. The details of such calculation are         well known to those skilled in the art.         Step 103 a. In this step a synchronization angle δ_(A) is         calculated as follows:     -   if there are phase-to-ground faults or phase-to-phase faults         then from the analyzing the total fault current flow FIG. 3,         FIG. 4, FIG. 5, FIG. 6, based on the assumption that faults         occurs in the line section between station A and series         capacitor and metal oxide varistor device SC&MOV the         synchronization angle δ_(A) in term e^(jδ) ^(A) from the formula         is determined:

$\begin{matrix} {\left\lbrack {\mathbb{e}}^{{j\delta}_{A}} \right\rbrack_{{{ph} - g},{{ph} - {ph}}} = \frac{{{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{I}}_{{SC}\; 2{\_ B}}} - {{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{I}}_{{SC}\; 1{\_ B}}}}{{{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{N}}_{A\; 1}} - {{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{N}}_{A\; 2}}}} & (1) \end{matrix}$ where: subscript ph-g denotes phase to ground faults: a-g, b-g, c-g, and subscript ph-ph denotes phase to phase faults a-b, b-c, c-a. a _(F2) ^(I-SET), a _(F1) ^(II-SET)—denote the coefficients dependent of fault type gathered in Table 1,

$\begin{matrix} {{{\underset{\_}{I}}_{{SC}\; 1{\_ B}} = {{{- \frac{1}{{\underset{\_}{Z}}_{c\; 1L}}}{{\sinh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{B\; 1}}} + {{\cosh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 1}}}},} & (2) \end{matrix}$ I _(SC1) _(—) _(B)—denotes the positive sequence current from the station B after analytical transfer towards series capacitor banks SC&MOV,

$\begin{matrix} {{{\underset{\_}{I}}_{{SC2}{\_ B}} = {{{- \frac{1}{{\underset{\_}{Z}}_{c\; 1L}}}{{\sinh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{B\; 2}}} + {{\cosh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 2}}}},} & (3) \end{matrix}$ I _(SC2) _(—) _(B)—denotes the negative sequence current from the station B after analytical transfer towards series capacitor banks SC&MOV,

$\begin{matrix} {{{\underset{\_}{N}}_{A\; 1} = {{{\cosh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)} \cdot {\underset{\_}{I}}_{A\; 1}} - {\frac{1}{{\underset{\_}{Z}}_{c\; 1L}}{{\sinh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)} \cdot {\underset{\_}{V}}_{A\; 1}}}}},} & (4) \\ {{{\underset{\_}{N}}_{A\; 2} = {{{\cosh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)} \cdot {\underset{\_}{I}}_{A\; 2}} - {\frac{1}{{\underset{\_}{Z}}_{c\; 1L}}{{\sinh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)} \cdot {\underset{\_}{V}}_{A\; 2}}}}},} & (5) \end{matrix}$ where: l—denotes the total line length, p_(SC)—denotes the distance (p.u.) per unit length, at which the compensation bank (SC&MOV) is installed far away from the station A, which is not presented on the drawings,

${\underset{\_}{Z}}_{c\; 1L} = {\sqrt{\frac{{\underset{\_}{Z}}_{1L}^{\prime}}{{\underset{\_}{Y}}_{1L}^{\prime}}}\mspace{14mu}\text{-}}$ denotes the characteristic impedance of the line for the positive sequence, γ _(1L)=√{square root over (Z′_(1L) Y′_(1L))}—denotes the propagation constant of the line for the positive sequence, Z′_(1L), Y′_(1L)—denote the positive sequence impedance and admittance of the line per unit length. It was considered for determining the negative-sequence quantities according to equations (3) and (5) that the positive- and negative-sequence line parameters are identical, as in reality. For both sequences the line parameters for the positive sequence (subscript: 1) are used. This is also applied in the following considerations. Analysis of boundary conditions for different types of faults shows that there is a certain degree of freedom when determining the share coefficients determining the relation between the symmetrical components of the total fault current when voltage drop across the fault resistance is estimated. Their selection depends on the adopted preference of the use of individual sequence components depending on the type of the fault. In the presented example of the invention embodiment, in order to ensure high precision of the fault location, voltage drop across the fault resistance is estimated using positive and negative sequence components.

There are two characteristic alternative sets (among the other possible) of the share coefficients for the phase-to-ground a-g, b-g, c-g, and phase-to-phase faults a-b, b-c, c-a, as gathered in the table 1.

TABLE 1 Phase-to-ground faults and phase-to-phase faults - two alternative sets of share coefficients. I-SET II-SET FAULT a _(F1) ^(I-SET) a _(F2) ^(I-SET) a _(F1) ^(II-SET) a _(F2) ^(II-SET) a-g 0 3 3 0 b-g 0 3a 3a ² 0 c-g 0 3a ² 3a 0 a-b 0 1 − a  1 − a ² 0 b-c 0  a − a ² a ² − a  0 c-a 0 a ² − 1  a − 1  0 a = exp(j2π/3) In contrast to the above fault types, the other situation is for the remaining fault types (phase-to-phase-to-ground and three phase symmetrical faults). This is so, since for these remaining faults there is no alternative sets of the share coefficients for the positive and the negative sequence—as gathered in the table 2.

TABLE 2 Set of share coefficients for phase-to-phase-to-ground faults and three phase faults. FAULT a _(F1) a _(F2) a-b-g, a-b-c,  1 − a ² 1 − a a-b-c-g b-c-g a ² − a   a − a ² c-a-g a − 1 a ² − 1  The type of the fault is denoted by symbols: a-g, b-g, c-g, a-b, b-c, c-a, a-b-g, a-b-c, a-b-c-g, b-c-g, c-a-g, where letters a, b, c denote individual phases, and letter g denotes earthing (ground), index 1 denotes the positive-sequence component, index 2—the negative component. If there are not phase-to-ground faults and not phase-to-phase faults, then the synchronization angle δ_(A) in term e^(jδ) ^(A) from the analysis of the current flow for the pre-fault conditions (FIG. 2) with using the distributed parameter line model is calculated from the formula:

$\begin{matrix} {{{{\mathbb{e}}^{{j\delta}_{A}} = {\frac{{\left( {1 + {{\underset{\_}{Z}}_{1B}{\underset{\_}{Y}}_{1B}}} \right)\left( {{{\underset{\_}{Y}}_{1B}{\underset{\_}{V}}_{B\; 1}^{pre}} - {\underset{\_}{I}}_{B\; 1}^{pre}} \right)} + {{\underset{\_}{Y}}_{1B}{\underset{\_}{V}}_{B\; 1}^{pre}}}{{\left( {1 + {{\underset{\_}{Z}}_{1A}{\underset{\_}{Y}}_{1A}}} \right)\left( {{\underset{\_}{I}}_{A\; 1}^{pre} - {{\underset{\_}{Y}}_{1A}{\underset{\_}{V}}_{A\; 1}^{pre}}} \right)} - {{\underset{\_}{Y}}_{1A}{\underset{\_}{V}}_{A\; 1}^{pre}}}\mspace{14mu}{where}\text{:}}}{\underset{\_}{Z}}_{1A} = {p_{SC}l{\underset{\_}{A}}_{1{sh\_ p}}{\underset{\_}{Z}}_{1L}^{\prime}}},{{\underset{\_}{Y}}_{1A} = {0.5p_{SC}l{\underset{\_}{A}}_{1{th\_ p}}{\underset{\_}{Y}}_{1L}^{\prime}}},{{\underset{\_}{Z}}_{1B} = {\left( {1 - p_{SC}} \right)l{\underset{\_}{A}}_{{1{sh\_}1} - p}{\underset{\_}{Z}}_{1L}^{\prime}}},{{\underset{\_}{Y}}_{1B} = {0.5\left( {1 - p_{SC}} \right)l{\underset{\_}{A}}_{{1{th\_}1} - p}{\underset{\_}{Y}}_{1L}^{\prime}}},{{\underset{\_}{A}}_{1{sh\_ p}} = \frac{\sinh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)}{{\underset{\_}{\gamma}}_{1L}p_{SC}\ell}},{{\underset{\_}{A}}_{1{th\_ p}} = \frac{\tanh\left( {0.5\;{\underset{\_}{\gamma}}_{1L}p_{sc}l} \right)}{0.5{\underset{\_}{\gamma}}_{1L}p_{SC}\ell}},{{\underset{\_}{A}}_{{1{sh\_}1} - p} = \frac{\sinh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)}{{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}\ell}},{{\underset{\_}{A}}_{{1{th\_}1} - p} = \frac{\tanh\left( {0.5\;{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{sc}} \right)}l} \right)}{0.5{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l}},} & (6) \end{matrix}$ V _(A1) ^(pre), V _(B1) ^(pre)—denote pre-fault positive sequence quantities (superscript ‘pre’ and subscript ‘1’) of voltage measured at the ends A and B, respectively, I _(A1) ^(pre), I _(B1) ^(pre)—denote pre-fault positive sequence quantities (superscript ‘pre’ and subscript ‘1’) of current measured of the faulted line at the ends A and B, respectively. Step 104 a. Having determined the synchronization angle δ_(A), and in consequence the synchronization operator e^(jδ) ^(A) , one multiplies the phasors of the signals from the station A by this operator. This assures common reference for the phasors of signals acquired asynchronously at both line ends. Taking into account the distributed parameter line model, the following generalised fault loop model is applied, covering different fault types, what is reflected with use of the appropriate weighting coefficients: a _(i) and share coefficients: a _(Fi), both dependent on a fault type, where i=0, 1, 2 and denote i-th symmetrical component, where the a ₁, a ₂, a ₀, a _(m0)—weighting coefficients gathered in the table 3 and a _(F1), a _(F2)—share coefficients gathered in the table 4.

$\begin{matrix} {{{\sum\limits_{i = 0}^{2}\;{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Ai}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Ai}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}}} \right)} \cdot {\mathbb{e}}^{{j\delta}\; A}}} - {R_{FA}{\sum\limits_{i = 1}^{2}\;\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LA}} \right)}p_{SC}l} \right)}}}} = 0} & (7) \end{matrix}$ for single line,

$\begin{matrix} {{{\sum\limits_{i = 0}^{2}\;{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Ai}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Ai}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}}} \right)} \cdot {\mathbb{e}}^{{j\delta}\; A}}} - {R_{FA}{\sum\limits_{i = 1}^{2}\;\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LA}} \right)}p_{SC}l} \right)}}} - {{\underset{\_}{a}}_{m\; 0}d_{LA}p_{SC}{\underset{\_}{Z}}_{0\; m}{{\underset{\_}{I}}_{{{Aparal}\_}0} \cdot {\mathbb{e}}^{{j\delta}_{A}}}}} = 0} & (8) \end{matrix}$ for double-circuit line where: d_(LA)—denotes the unknown hypothetical distance to fault in [p.u.], expressed in relation to the length (p_(SC)l) of the line section between the station A and the series capacitor bank, R_(FA)—denotes the unknown fault resistance,

${{\underset{\_}{M}}_{1} = {{{\cosh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)}{\underset{\_}{I}}_{A\; 1}{\mathbb{e}}^{{j\delta}_{A}}} - \frac{{\sinh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)}{\underset{\_}{V}}_{A\; 1}{\mathbb{e}}^{{j\delta}_{A}}}{{\underset{\_}{Z}}_{c\; 1L}} - \frac{{\sinh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{B\; 1}}{{\underset{\_}{Z}}_{c\; 1L}} + {{\cosh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 1}}}},{{\underset{\_}{M}}_{2} = {{{\cosh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)}{\underset{\_}{I}}_{A\; 2}{\mathbb{e}}^{{j\delta}_{A}}} - \frac{{\sinh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)}{\underset{\_}{V}}_{A\; 2}{\mathbb{e}}^{{j\delta}_{A}}}{{\underset{\_}{Z}}_{c\; 1L}} - \frac{{\sinh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{B\; 2}}{{\underset{\_}{Z}}_{c\; 1L}} + {{\cosh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 2}}}},$ a ₁, a ₂, a ₀, a _(m0)—denote the weighting coefficients dependent on fault type gathered in the table 3, a _(F1), a _(F2)—denote the share coefficients dependent on fault type gathered in the table 4, Z _(0m)—mutual coupling impedance for zero sequence, I _(Aparal) _(—) ₀—denotes zero sequence current from the station A measured in healthy parallel line

TABLE 3 Weighting Coefficients for Composing Fault loop Voltage And Current FAULT a ₁ a ₂ a ₀ a _(m0) a-g 1 1 1 1 b-g −0.5 − j0.5{square root over (3)}  0.5 + j0.5{square root over (3)} 1 1 c-g  0.5 + j0.5{square root over (3)} −0.5 − j0.5{square root over (3)} 1 1 a-b, a-b-g  1.5 + j0.5{square root over (3)}  1.5 − j0.5{square root over (3)} 0 0 a-b-c, a-b-c-g b-c, b-c-g  −j{square root over (3)} j{square root over (3)} 0 0 c-a, c-a-g −1.5 + j0.5{square root over (3)} −1.5 − j0.5{square root over (3)} 0 0 a _(0m) - used in case of double circuit lines for mutual coupling compensation

TABLE 4 Share Coefficients Used For calculation Of Total Fault Current FAULT a _(F1) a _(F2) a _(F0) a-g 0 3 0 b-g 0  1.5 + j1.5{square root over (3)} 0 c-g 0 −1.5 − j1.5{square root over (3)} 0 a-b 0  1.5 − j0.5{square root over (3)} 0 b-c 0 j{square root over (3)} 0 c-a 0 −1.5 − j0.5{square root over (3)} 0 a-b-g 1.5 + j0.5{square root over (3)}  1.5 − j0.5{square root over (3)} 0 b-c-g −j{square root over (3)} j{square root over (3)} 0 c-a-g −1.5 + j0.5{square root over (3)}  −1.5 − j0.5{square root over (3)} 0 a-b-c 1.5 + j0.5{square root over (3)}   1.5 − j0.5{square root over (3)}*⁾ 0 a-b-c-g j = {square root over (−1)}, *⁾ a _(F2) ≠ 0, however the negative sequence components are not present in the signals Resolving equation (7) or (8) into the real/imaginary parts allows to solve it for the unknowns: d_(LA) and R_(FA). As for example, the iterative Newton-Raphson method can be applied for that. This method requires setting the starting values for the sought values. These values can be set to: 0.5 [p.u.] for d_(LA), and 0 for R_(FA). After completing the solution of (7) or (8), the distance to fault d_(A) [p.u.], which is related to the whole line length l, between the stations A and B, is calculated as: d _(A) =p _(sc) d _(LA)  (9) where d_(LA), is the result from the last iterative solution of the equation (7) or (8). Step 105 a. In this step the equivalent impedance of compensating bank at fault stage is calculated with using the post fault values of measured voltage and currents.

$\begin{matrix} {{{\underset{\_}{Z}}_{{SC}\; 1{\_{ph}}}^{SUB\_ A} = {\frac{\Delta\;{\underset{\_}{V}}_{SC\_ ph}^{SUB\_ A}}{{\underset{\_}{I}}_{SC\_ ph}^{SUB\_ A}}\mspace{14mu}{where}\text{:}}}{{{\Delta{\underset{\_}{V}}_{SC\_ ph}^{SUB\_ A}} = {{{\underset{\_}{a}}_{1\_}\Delta{\underset{\_}{V}}_{{SC\_}\; 1}^{SUB\_ A}} + {{\underset{\_}{a}}_{2\_}\Delta{\underset{\_}{V}}_{{SC\_}\; 2}^{SUB\_ A}} + {{\underset{\_}{a}}_{0\_}\Delta{\underset{\_}{V}}_{{SC\_}\; 0}^{SUB\_ A}}}},{{\underset{\_}{I}}_{SC\_ ph}^{SUB\_ A} = {{{\underset{\_}{a}}_{1\_}\Delta{\underset{\_}{I}}_{{SC\_}\; 1}^{SUB\_ A}} + {{\underset{\_}{a}}_{2\_}\Delta{\underset{\_}{I}}_{{SC\_}\; 2}^{SUB\_ A}} + {{\underset{\_}{a}}_{0\_}{\underset{\_}{I}}_{{SC\_}\; 0}^{SUB\_ A}}}},}} & (10) \end{matrix}$ ΔV _(SC) _(—) _(ph) ^(SUB) ^(—) ^(A), I _(SC) _(—) _(ph) ^(SUB) ^(—) ^(A)—denote the phase quantities of voltage and currents calculated from i-th symmetrical components of voltage drop and i-th symmetrical component of current entering SC&MOV bank, with neglecting the shunt capacitances of the line:

-   -   in case of single phase to ground fault: the faulted phase     -   in case of inter-phase fault: any of the faulted phases,         ΔV _(SC) _(—i) ^(SUB) ^(—) ^(A)=V _(Bi)−V _(Ai)e^(jδ) ^(A)         −Z_(iL)I_(Bi)+d_(A)(Z _(iL) I _(Ai)e^(jδ) ^(A)         +Z_(iL)I_(Bi))—denotes the i-th symmetrical component of voltage         drop across the SCs&MOVs,         I _(SC) _(—) _(i) ^(SUB) ^(—) ^(A)=I _(Bi)—denotes i-th         symmetrical component of current entering SC&MOV bank, with         neglecting the shunt capacitances of the line,         a ₁ _(—) , a ₂ _(—) , a ₀ _(—) —complex coefficients dependent         on fault type which are selected in such a way that in case of:     -   single phase to ground fault—the faulted phase, as for example         for a-g fault the equivalent impedance of SC&MOV in phase ‘a’ is         determined,     -   inter-phase fault—any of the faulted phases, as for example for         a-b fault the equivalent impedance of SC&MOV in phase ‘a’ is         determined.         Step 106 a. At this stage the equivalent impedance of         compensation bank at pre-fault stage is calculated:

$\begin{matrix} {\mspace{79mu}{{{\underset{\_}{Z}}_{{SC}\; 1{\_ pre}} = {\frac{\Delta{\underset{\_}{V}}_{{SC}\; 1}^{pre}}{{\underset{\_}{I}}_{{SC}\; 1{\_ B}}^{pre}}\mspace{14mu}{where}\text{:}}}\mspace{79mu}{{{\Delta{\underset{\_}{V}}_{{SC}\; 1}^{pre}} = {{\underset{\_}{V}}_{{SC}\; 1{\_ B}}^{pre} - {{\underset{\_}{V}}_{{SC}\; 1{\_ A}}^{pre}{\mathbb{e}}^{{j\delta}_{A}}}}},\mspace{79mu}{{\underset{\_}{V}}_{{SC}\; 1{\_ A}}^{pre} = {{{\cosh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)} \cdot {\underset{\_}{V}}_{A\; 1}^{pre}} - {{\underset{\_}{Z}}_{c\; 1L}{{\sinh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)} \cdot {\underset{\_}{I}}_{A\; 1}^{pre}}}}},{{\underset{\_}{V}}_{{SC}\; 1{\_ B}}^{pre} = {{{\cosh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{A\; 1}^{pre}} - {{\underset{\_}{Z}}_{c\; 1L}{{\sinh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 1}^{pre}}}}},{{\underset{\_}{I}}_{{SC}\; 1{\_ B}}^{pre} = {{{- \frac{1}{{\underset{\_}{Z}}_{c\; 1L}}}{{\sinh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{B\; 1}^{pre}}} + {{\cosh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {{\underset{\_}{I}}_{B\; 1}^{pre}.}}}}}}} & (11) \end{matrix}$ Step 103 b. In this step the synchronization angle is calculated as follows: if there is phase-to-ground faults or phase-to-phase faults then from the analyzing the total fault current flow (FIG. 7, FIG. 8, FIG. 9, FIG. 10), based on the assumption that faults occurs in the line section between station B and series capacitor and metal oxide varistor device SC&MOV determining the synchronization angle δ_(B) in term e^(jδ) ^(B) from the formula:

$\begin{matrix} {\left\lbrack {\mathbb{e}}^{{j\delta}_{B}} \right\rbrack_{{{ph} - g},{{ph} - {ph}}} = \frac{{{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{I}}_{{SC}\; 2{\_ A}}} - {{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{I}}_{{SC}\; 1{\_ A}}}}{{{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{N}}_{B\; 1}} - {{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{N}}_{B\; 2}}}} & (12) \end{matrix}$ where: a _(F2) ^(I-SET), a _(F1) ^(II-SET)—denote the coefficients dependent of fault type gathered in table 1, F2

$\begin{matrix} {{{\underset{\_}{I}}_{{SC}\; 1{\_ A}} = {{{- \frac{1}{{\underset{\_}{Z}}_{c\; 1L}}}{{\sinh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)} \cdot {\underset{\_}{V}}_{A\; 1}}} + {{\cosh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)} \cdot {\underset{\_}{I}}_{A\; 1}}}},} & (13) \end{matrix}$

I _(SC1) _(—) _(A)—denotes positive sequence current from the station A after analytical transfer towards series capacitor banks SC&MOV,

$\begin{matrix} {{{\underset{\_}{I}}_{{SC}\; 2{\_ A}} = {{{- \frac{1}{{\underset{\_}{Z}}_{c\; 1L}}}{{\sinh\left( {{\underset{\_}{\gamma}}_{1\; L}p_{SC}l} \right)} \cdot {\underset{\_}{V}}_{A\; 2}}} + {{\cosh\left( {{\underset{\_}{\gamma}}_{1L}p_{SC}l} \right)} \cdot {\underset{\_}{I}}_{A\; 2}}}},} & (14) \end{matrix}$ I _(SC2) _(—) _(A)—denotes negative sequence current from the station A after analytical transfer towards series capacitor banks SC&MOV,

$\begin{matrix} {{{\underset{\_}{N}}_{B\; 1} = {{{\cosh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 1}} - {\frac{1}{{\underset{\_}{Z}}_{c\; 1L}}{{\sinh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{B\; 1}}}}},} & (15) \\ {{{\underset{\_}{N}}_{B\; 2} = {{{\cosh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 2}} - {\frac{1}{{\underset{\_}{Z}}_{c\; 1L}}{{\sinh\left( {{{\underset{\_}{\gamma}}_{1L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{B\; 2}}}}},} & (16) \end{matrix}$ where: l—denotes total line length, p_(SC)—denotes distance (p.u.) at which the compensation bank (SCs&MOVs) is installed far away from the station A,

${\underset{\_}{Z}}_{c\; 1L} = \sqrt{\frac{{\underset{\_}{Z}}_{1L}^{\prime}}{{\underset{\_}{Y}}_{1L}^{\prime}}}$ denotes characteristic impedance of the line for the positive sequence, γ _(1L)=√{square root over (Z′_(1L) Y′_(1L))}—denotes propagation constant of the line for the positive sequence, Z′_(1L) Y′_(1L)—denote positive sequence impedance and admittance of the line per unit length. It was considered for determining the negative-sequence quantities according to equations (14) and (16) that the positive- and negative-sequence line parameters are identical, as in reality. For both sequences the line parameters for the positive sequence (subscript: 1) are used. This is also applied in the following considerations:

if there are not phase-to-ground faults and not phase-to-phase faults then determining the synchronization angle δ_(B) for phase-to-ground faults and phase-to-phase faults in term e^(jδ) ^(B) from the formula (6).

Step 104 b. Having determined the synchronization angle δ_(B), and in consequence the synchronization operator e^(jδ) ^(B) , one multiplies the phasors of the signals from the end B by this operator. This assures common reference for the phasors of signals acquired asynchronously at both line ends.

Taking into account the distributed parameter line model, the following generalised fault loop model, covering different fault types, what is reflected with use of the appropriate weighting coefficients: a _(i) and share coefficients: a _(Fi), both dependent on a fault type is applied, where the a ₁, a ₂, a ₀, a _(m0)—weighting coefficients gathered in the table 3 and a _(F1), a _(F2)—share coefficients gathered in the table 4.

$\begin{matrix} {{{\sum\limits_{i = 0}^{2}\;{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Bi}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Bi}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}}} \right)} \cdot {\mathbb{e}}^{{j\delta}_{B}}}} - {R_{FB}{\sum\limits_{i = 1}^{2}\;\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LB}} \right)}\left( {1 - p_{SC}} \right)l} \right)}}}} = 0} & (17) \end{matrix}$ for single line,

$\begin{matrix} {{{\sum\limits_{i = 0}^{2}\;{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Bi}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Bi}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}}} \right)} \cdot {\mathbb{e}}^{{j\delta}_{B}}}} - {R_{FB}{\sum\limits_{i = 1}^{2}\;\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LB}} \right)}\left( {1 - p_{SC}} \right)l} \right)}}} - {{\underset{\_}{a}}_{m\; 0}{d_{LB}\left( {1 - p_{SC}} \right)}{\underset{\_}{Z}}_{0\; m}{{\underset{\_}{I}}_{{Bparal\_}0} \cdot {\mathbb{e}}^{{j\delta}_{B}}}}} = 0} & (18) \end{matrix}$ for double-circuit line where: d_(LB)—denotes the unknown hiphotetical distance to fault in [p.u.], expressed in relation to the length ((1−P_(SC))l) of the line section between the station B and the series capacitor bank, R_(FB)—denotes the unknown fault resistance,

${{\underset{\_}{M}}_{1} = {{{\cosh\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)}{\underset{\_}{I}}_{B\; 1}{\mathbb{e}}^{{j\delta}_{B}}} - \frac{{\sinh\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)}{\underset{\_}{V}}_{B\; 1}{\mathbb{e}}^{{j\delta}_{B}}}{{\underset{\_}{Z}}_{c\; 1\; L}} - \frac{{\sinh\left( {{\underset{\_}{\gamma}}_{1\; L}p_{SC}l} \right)} \cdot {\underset{\_}{V}}_{A\; 1}}{{\underset{\_}{Z}}_{c\; 1\; L}} + {{\cosh\left( {{\underset{\_}{\gamma}}_{1\; L}p_{SC}l} \right)} \cdot {\underset{\_}{I}}_{A\; 1}}}},{{\underset{\_}{M}}_{2} = {{{\cosh\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)}{\underset{\_}{I}}_{B\; 2}{\mathbb{e}}^{{j\delta}_{B}}} - \frac{{\sinh\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)}{\underset{\_}{V}}_{B\; 2}{\mathbb{e}}^{{j\delta}_{B}}}{{\underset{\_}{Z}}_{c\; 1\; L}} - \frac{{\sinh\left( {{\underset{\_}{\gamma}}_{1\; L}p_{SC}l} \right)} \cdot {\underset{\_}{V}}_{A\; 2}}{{\underset{\_}{Z}}_{c\; 1\; L}} + {{\cosh\left( {{\underset{\_}{\gamma}}_{1\; L}p_{SC}l} \right)} \cdot {\underset{\_}{I}}_{A\; 2}}}},$ a ₁, a ₂, a ₀, a _(m0)—denote weighting coefficients dependent on fault type gathered in table 3, a _(F1), a _(F2)—denote share coefficients dependent on fault type, gathered in table 4, Z _(0m)—mutual coupling impedance for zero sequence, I _(Bparal) _(—) ₀—denotes zero sequence current from the station B measured in parallel line Resolving equation (17) or (18) into the real/imaginary parts allows to solve it for the unknowns: d_(LB), R_(FB). As for example, the iterative Newton-Raphson method can be applied for that. This method requires setting the starting values for the sought values. These values can be set to: 0.5 [p.u.] for d_(LB) and 0 for R_(FB). After completing the solution of (17) or (18), the distance to fault d_(B) [p.u.], which is related to the whole line length l, between the stations A and B, is calculated as: d _(B)=(1−p _(SC))d _(LB)  (19) where d_(LB) is the result from the last iterative solution of (17) or (18). Step 105B. The equivalent impedance of compensating bank at fault stage is calculated with using the post fault values of measured voltage and currents:

$\begin{matrix} {{{\underset{\_}{Z}}_{{SC}\; 1{\_ ph}}^{SUB\_ B} = \frac{\Delta\;{\underset{\_}{V}}_{{SC}{\_ ph}}^{SUB\_ B}}{{\underset{\_}{I}}_{{SC}{\_ ph}}^{SUB\_ B}}}{{where}\text{:}}{{{\Delta{\underset{\_}{V}}_{{SC}{\_ ph}}^{SUB\_ B}} = {{{\underset{\_}{a}}_{1\_}\Delta{\underset{\_}{V}}_{{{SC}\_}\; 1}^{SUB\_ B}} + {{\underset{\_}{a}}_{2\_}\Delta{\underset{\_}{V}}_{{{SC}\_}\; 2}^{SUB\_ B}} + {{\underset{\_}{a}}_{0\_}\Delta{\underset{\_}{V}}_{{{SC}\_}\; 0}^{SUB\_ B}}}},{{\underset{\_}{I}}_{{SC}{\_ ph}}^{SUB\_ B} = {{{\underset{\_}{a}}_{1\_}\Delta{\underset{\_}{I}}_{{{SC}\_}\; 1}^{SUB\_ B}} + {{\underset{\_}{a}}_{2\_}\Delta{\underset{\_}{I}}_{{{SC}\_}\; 2}^{SUB\_ B}} + {{\underset{\_}{a}}_{0\_}{\underset{\_}{I}}_{{{SC}\_}\; 0}^{SUB\_ B}}}},}} & (20) \end{matrix}$ ΔV _(SC) _(—) _(ph) ^(SUB) ^(—) ^(A), I _(SC) _(—) _(ph) ^(SUB) ^(—) ^(A)—denote phase quantities of voltage and currents calculated from i-th symmetrical components of voltage drop and i-th symmetrical component of current entering SC&MOV bank, with neglecting the shunt capacitances of the line

-   -   in case of single phase to ground fault: the faulted phase,     -   in case of inter-phase fault: any of the faulted phases,         a ₁ _(—) , a ₂ _(—) , a ₀ _(—) —denote complex coefficients         dependent on fault type, which are selected in such a way that         in case of:     -   single phase to ground fault: the faulted phase, as for example         for a-g fault the equivalent impedance of SC&MOV in phase ‘a’ is         determined,     -   inter-phase fault: any of the faulted phases, as for example for         a-b fault the equivalent impedance of SC&MOV in phase ‘a’ is         determined,         ΔV _(SC) _(—i) ^(SUB) ^(—) ^(B)=V _(Ai)−V _(Bi)e^(jδ) ^(B)         −Z_(iL)I_(Bi)+d_(B)(Z _(iL) I _(Ai)e^(jδ) ^(A)         +Z_(iL)I_(Bi)e^(jδ) ^(B) )—denotes the voltage drop across the         SCs&MOVs for the i-th symmetrical component,         I _(SC) _(—) _(i) ^(SUB) ^(—) ^(B)=I _(Ai)—denote the i-th         symmetrical component of current entering SC&MOV bank, with         neglecting the shunt capacitances of the line.         Step 106 b. At this stage the equivalent impedance of         compensation bank at pre-fault stage is calculated as in step         106 a         Step 107. At this step the selection of final results is         generated.         First, it is check whether the particular subroutine I and         subroutine II yields the distance to fault within the section         range and whether the fault resistance is of positive value. If         this is not, so (if the determined distance to fault lays         outside the section range and/or the fault resistance is         negative) then such subroutine is rejected.         Next, the features “real” and:“imag” of the estimated equivalent         circuit of compensation banks SCs&MOVS are checked.         If it is satisfied the following dependence:         real(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A))>0,         imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A))<0         and imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A))>imag(Z _(SC1)         _(—) _(pre))         that denotes, that the fault has occurred between the station A         and the SC&MOV bank, and the results:         d_(A)—distance to fault (p.u.), counted from the station A up to         fault point F (FIG. 1),         R_(FA)—fault resistance,         are the valid results, i.e. consistent with the real fault case.         If opposite dependences are satisfying, that means:     -   real(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B))>0,     -   imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B))<0     -   and imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B))>imag(Z _(SC1)         _(—) _(pre))     -   that denotes that the fault has occurred between the station B         and the SC&MOV bank, and the results:     -   d_(B)—distance to fault (p.u.), counted from the station B up to         fault point F (FIG. 1),     -   R_(FB)—fault resistance,     -   are the valid results, i.e. consistent with the real fault case.         It is observed that selection of d_(A) or d_(B) is independent         of the status and/or value of series capacitance. 

The invention claimed is:
 1. A computer program product comprising non-transitory computer-readable program code which when executed on a computing device carries out the steps of detecting a location of a fault in one or more power transmission lines that extend between a station A and a station B and have a series capacitor and metal oxide varistor device (SC&MOV) connected therein, using unsynchronized voltage and current measurements in the stations A and B, before and after occurrence of the fault, performing a subroutine I for determining if the fault occurred in a line section (LA) between the station A and the SC&MOV and determining a distance (d_(A)) to the fault from station A, performing a subroutine II for determining if the fault occurred in a line section (LB) between the station B and the SC&MOV and determining a distance (d_(B)) to the fault from station B wherein subroutine I comprises the following steps: (A1.) determining a synchronization angle (δA) in term (e^(jδ) ^(A) ) for a known type of fault which is phase-to-ground faults or phase-to-phase faults, from the formula: $\left\lbrack {\mathbb{e}}^{{j\delta}_{A}} \right\rbrack_{{{ph} - g},{{ph} - {ph}}} = \frac{{{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{I}}_{{SC}\; 2{\_ B}}} - {{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{I}}_{{SC}\; 1{\_ B}}}}{{{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{N}}_{A\; 1}} - {{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{N}}_{A\; 2}}}$ where: subscript ph-g denotes phase to ground faults: a-g, b-g, c-g, and subscript ph-ph denotes phase to phase faults a-b, b-c, c-a, a _(F2) ^(I-SET), a _(F1) ^(II-SET)—denote coefficients dependent on fault type, I _(SC1) _(—) _(B)—denotes the positive sequence current from the station B after analytical transfer towards the SC&MOV, I _(SC2) _(—) _(B)—denotes the negative sequence current from the station B after analytical transfer towards the SC&MOV, N _(A1), N _(A2)—denote coefficients dependent on positive and negative sequence currents and voltage taken at station A respectively as well as parameters of the one or more power transmission lines, or for faults not being the phase-to-ground faults or phase-to-phase faults from the formula: ${\mathbb{e}}^{{j\delta}_{A}} = \frac{{\left( {1 + {{\underset{\_}{Z}}_{1\; B}{\underset{\_}{Y}}_{1\; B}}} \right)\left( {{{\underset{\_}{Y}}_{1\; B}{\underset{\_}{V}}_{B\; 1}^{pre}} - {\underset{\_}{I}}_{B\; 1}^{pre}} \right)} + {{\underset{\_}{Y}}_{1\; B}{\underset{\_}{V}}_{B\; 1}^{pre}}}{{\left( {1 + {{\underset{\_}{Z}}_{1\; A}{\underset{\_}{Y}}_{1\; A}}} \right)\left( {{\underset{\_}{I}}_{A\; 1}^{pre} - {{\underset{\_}{Y}}_{1\; A}{\underset{\_}{V}}_{A\; 1}^{pre}}} \right)} - {{\underset{\_}{Y}}_{1\; A}{\underset{\_}{V}}_{A\; 1}^{pre}}}$ where: V _(A1) ^(pre), V _(B1) ^(pre)—denote pre-fault positive sequence quantities of voltage measured at the ends A and B, respectively, I _(A1) ^(pre), I _(B1) ^(pre)—denote pre-fault positive sequence quantities of current measured of the faulted one or more power transmission lines at the ends A and B, respectively, Z _(1A), Y _(1A)—denote impedance and admittance, respectively, of line section LA, Z_(1B), Y_(1B)—denote impedance and admittance, respectively, of line section LB, (A2.) then taking into account a distributed parameter line model and solving the formula: d _(A) =p _(SC) d _(LA), where: p_(sc)—denotes a distance (p.u.) per unit of total line length l at which the SC&MOV is installed in the transmission line from the station A, and d_(LA)—denotes a hypothetical distance to fault, expressed in relation to length p_(SC)l and which value is received using known iterative method by solving the fault loop formula: ${{\sum\limits_{i = 0}^{2}\;{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Ai}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Ai}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}}} \right)} \cdot {\mathbb{e}}^{{j\delta}_{A}}}} - {R_{FA}{\sum\limits_{i = 1}^{2}\;\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LA}} \right)}p_{SC}l} \right)}}}} = 0$ for single line, in which: P_(SC)l—denotes a length (km) from station A at which SC&MOV is installed in the transmission line, V _(Ai)—zero, positive, and negative voltage sequences measured at the end A, and wherein i=0, 1, 2, denoting the i-th symmetrical component, Z _(ciL)—zero, positive, and negative sequence line impedances, and wherein i=0, 1, 2, denoting the i-th symmetrical component, I _(Ai)—zero, positive, and negative sequence currents measured at the end A, and wherein i=0, 1, 2, denoting the i-th symmetrical component, γ _(iL)—zero, positive, and negative sequences of line admittance, and wherein i=0, 1, 2, denoting the i-th symmetrical component; ${{\sum\limits_{i = 0}^{2}\;{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Ai}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Ai}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}d_{LA}p_{SC}l} \right)}}} \right)} \cdot {\mathbb{e}}^{{j\delta}_{A}}}} - {R_{FA}{\sum\limits_{i = 1}^{2}\;\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}\left( {1 - {d_{LA}p_{SC}l}} \right)} \right.}}} - {{\underset{\_}{a}}_{m\; 0}d_{LA}p_{SC}{\underset{\_}{Z}}_{0\; m}{{\underset{\_}{I}}_{{A{paral\_}}0} \cdot {\mathbb{e}}^{{j\delta}_{A}}}}} = 0$ for double-circuit line, in which: R_(FA)—denotes an unknown fault resistance, which a value is received using known iterative method by solving the fault loop formula above, a ₁, a ₂, a ₀, a _(m0)—denotes weighting coefficients dependent on fault type, M _(i)—denote numerical components for positive and negative sequences, a _(F0), a _(F1), a _(F2)—denotes the share coefficients dependent on fault type, Z _(0m)—mutual coupling impedance for zero sequence, I _(Aparal) _(—) ₀—denotes zero sequence current from the station A measured in a healthy parallel line, V _(Ai)—zero, positive, and negative voltage sequences measured at the end A, and wherein i=0, 1, 2, denoting the i-th symmetrical component, Z _(ciL)—zero, positive, and negative sequence line impedances, and wherein i=0, 1, 2, denoting the i-th symmetrical component, I _(Ai)—zero, positive, and negative sequence currents measured at the end A, and wherein i=0, 1, 2, denoting the i-th symmetrical component, γ _(iL)—zero, positive, and negative sequences of line admittance, and wherein i=0, 1, 2, denoting the i-th symmetrical component, (A3.) then an equivalent impedance of the SC&MOV at fault stage (Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A)) is calculated from the equation: ${\underset{\_}{Z}}_{{SC}\; 1{\_ ph}}^{SUB\_ A} = \frac{\Delta\;{\underset{\_}{V}}_{SC\_ ph}^{SUB\_ A}}{{\underset{\_}{I}}_{SC\_ ph}^{SUB\_ A}}$ where: ${{\Delta\;{\underset{\_}{V}}_{SC\_ ph}^{SUB\_ A}} = {{{\underset{\_}{a}}_{1}\Delta\;{\underset{\_}{V}}_{{SC\_}1}^{SUB\_ A}} + {{\underset{\_}{a}}_{2}\Delta\;{\underset{\_}{V}}_{{SC\_}2}^{SUB\_ A}} + {{\underset{\_}{a}}_{0}\Delta\;{\underset{\_}{V}}_{{SC\_}0}^{SUB\_ A}}}},{{\underset{\_}{I}}_{SC\_ ph}^{SUB\_ A} = {{{\underset{\_}{a}}_{1}\Delta\;{\underset{\_}{I}}_{{SC\_}1}^{SUB\_ A}} + {{\underset{\_}{a}}_{2}\Delta\;{\underset{\_}{I}}_{{SC\_}2}^{SUB\_ A}} + {{\underset{\_}{a}}_{0}{\underset{\_}{I}}_{{SC\_}0}^{SUB\_ A}}}},$ ΔV _(SC) _(—) _(ph) ^(SUB) ^(—) ^(A) and I _(SC) _(—) _(ph) ^(SUB) ^(—) ^(A) denote the phase quantities of voltage and currents calculated from the i-th symmetrical components of voltage drop and i-th symmetrical component of current entering the SC&MOV bank while neglecting the shunt capacitances of the line, ΔV _(SC) _(—i) ^(SUB) ^(—) ^(A)=V _(Bi)−V _(Ai)e^(jδ) ^(A) −Z_(iL)I_(Bi)+d_(A)(Z _(iL) I _(Ai)e^(jδ) ^(A) +Z_(iL)I_(Bi)) denotes the i-th symmetrical component of voltage drop across the SCs&MOVs, I _(SC) _(—) _(i) ^(SUB) ^(—) ^(A)=I _(Bi) denotes i-th symmetrical component of current entering the SC&MOV bank, while neglecting the shunt capacitances of the line, a ₁, a ₂, a ₀ are complex coefficients that depend on fault type and are selected by case type: in a single-phase to ground fault, the faulted phase, for example, for a phase ‘a’ to ground fault the equivalent impedance of SC&MOV in phase ‘a’ is determined and in an inter-phase fault, any of the faulted phases, for example, for a phase ‘a’ to phase ‘b’ fault, the equivalent impedance of SC&MOV in phase ‘a’ is determined, next, an-equivalent impedance of SC&MOV at pre-fault (Z_(SC1) _(—) _(pre)) is calculated from the equation: $\mspace{79mu}{{\underset{\_}{Z}}_{{SC}\; 1{\_ pre}} = \frac{\Delta\;{\underset{\_}{V}}_{{SC}\; 1}^{pre}}{{\underset{\_}{I}}_{{SC}\; 1{\_ B}}^{pre}}}$      where: $\mspace{79mu}{{{\Delta\;{\underset{\_}{V}}_{{SC}\; 1}^{pre}} = {{\underset{\_}{V}}_{{SC}\; 1{\_ B}}^{pre} - {{\underset{\_}{V}}_{{SC}\; 1{\_ A}}^{pre}{\mathbb{e}}^{j\;\delta\; A}}}},\mspace{79mu}{{\underset{\_}{V}}_{{SC}\; 1{\_ A}}^{pre} = {{\cos\;{{h\left( {{\underset{\_}{\gamma}}_{1\; L}p_{SC}l} \right)} \cdot {\underset{\_}{V}}_{A\; 1}^{pre}}} - {{\underset{\_}{Z}}_{c\; 1L}\sin\;{{h\left( {{\underset{\_}{\gamma}}_{1\; L}p_{SC}l} \right)} \cdot {\underset{\_}{I}}_{A\; 1}^{pre}}}}},{{\underset{\_}{V}}_{{SC}\; 1{\_ B}}^{pre} = {{\cos\; h{\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right) \cdot {\underset{\_}{V}}_{B\; 1}^{pre}}} - {{\underset{\_}{Z}}_{c\; 1L}\sin\;{{h\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 1}^{pre}}}}},{{\underset{\_}{I}}_{{SC}\; 1{\_ B}}^{pre} = {{{- \frac{1}{{\underset{\_}{Z}}_{c\; 1\; L}}}\sin\;{{h\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{B\; 1}^{pre}}} + {\cos\;{{h\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 1}^{pre}}}}},}$ in which: Z_(c1L) denotes the characteristic impedance of the line for the positive sequence, and γ _(1L) denotes the propagation constant of the line for the positive sequence, in order to determine the final result (d_(A)) which is chosen in the following steps: first, check whether the fault resistances (R_(FA)) is of positive value and if not, then subroutine I is rejected, next, features “real” and “imag” of the estimated equivalent circuit, (Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A)), of the SC&MOV is checked for the following dependency: real(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A))>0, imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A))<0, and imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(A))>imag(Z _(SC1) _(—) _(pre)) and if the above dependency is met, then it is determined that the fault has occurred between the station A and the SC&MOV and the result (d_(A)) indicates the distance to fault, and the information about the real fault case is supplied to fault locator FL, wherein subroutine II comprises the following steps: (B1.) determining a synchronization angle (δ_(B)) in term (e^(jδ) ^(B) ) for a known type of fault which is phase-to-ground faults or phase-to-phase faults from the formula: $\left\lbrack {\mathbb{e}}^{{j\delta}_{B}} \right\rbrack_{{{ph} - g},{{ph} - {ph}}} = \frac{{{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{I}}_{{SC}\; 2{\_ A}}} - {{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{I}}_{{SC}\; 1{\_ A}}}}{{{\underset{\_}{a}}_{F\; 1}^{{II} - {SET}}{\underset{\_}{N}}_{B\; 1}} - {{\underset{\_}{a}}_{F\; 2}^{I - {SET}}{\underset{\_}{N}}_{B\; 2}}}$ where: a _(F2) ^(I-SET), a _(F1) ^(II-SET)—denote coefficients dependent of fault type, I _(SC1) _(—) _(A)—denotes positive sequence current from the station A after analytical transfer towards the SC&MOV, I _(SC2) _(—) _(A)—denotes negative sequence current from the station A after analytical transfer towards the SC&MOV, N _(B1), N _(B2)—denote coefficients dependent on positive and negative sequence currents and voltage taken at station B respectively as well as parameters of the one or more power transmission lines, and for faults not being the phase-to-ground faults or phase-to-phase faults from the formula: ${\mathbb{e}}^{- {j\delta}_{B}} = \frac{{\left( {1 + {{\underset{\_}{Z}}_{1\; B}{\underset{\_}{Y}}_{1\; B}}} \right)\left( {{{\underset{\_}{Y}}_{1\; B}{\underset{\_}{V}}_{B\; 1}^{pre}} - {\underset{\_}{I}}_{B\; 1}^{pre}} \right)} + {{\underset{\_}{Y}}_{1\; B}{\underset{\_}{V}}_{B\; 1}^{pre}}}{{\left( {1 + {{\underset{\_}{Z}}_{1\; A}{\underset{\_}{Y}}_{1\; A}}} \right)\left( {{\underset{\_}{I}}_{A\; 1}^{pre} - {{\underset{\_}{Y}}_{1\; A}{\underset{\_}{V}}_{A\; 1}^{pre}}} \right)} - {{\underset{\_}{Y}}_{1\; A}{\underset{\_}{V}}_{A\; 1}^{pre}}}$ where: V _(A1) ^(pre), V _(B1) ^(pre)—denote pre-fault positive sequence quantities of voltage measured at the ends A and B, respectively, I _(A1) ^(pre), I _(B1) ^(pre)—denote pre-fault positive sequence quantities of current measured of the faulted one or more power transmission lines at the ends A and B, respectively, Z _(1A), Y _(1A)—denote impedance and admittance, respectively, of line section LA, Z _(1B), Y _(1B)—denote impedance and admittance, respectively, of line section LB, (B2.) then taking into account the distributed parameter line model and solving the formula: d _(B)=(1−p _(SC))d _(LB) where: (1−p_(SC))—denotes a distance (p.u.) per unit of total line length l, at which the SC&MOV is installed in the transmission line from the station B, d_(LB)—denotes a hypothetical distance to fault, expressed in relation to length (1−p_(SC)l) and which value is received using known iterative method by solving the fault loop formula: ${{\sum\limits_{i = 0}^{2}\;{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Bi}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Bi}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}}} \right)} \cdot {\mathbb{e}}^{{j\delta}_{B}}}} - {R_{FB}{\sum\limits_{i = 1}^{2}\;\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LB}} \right)}\left( {1 - p_{SC}} \right)l} \right)}}}} = 0$ for single line, ${{\sum\limits_{i = 0}^{2}\;{{{\underset{\_}{a}}_{i}\left( {{{\underset{\_}{V}}_{Bi}{\cosh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}} - {{\underset{\_}{Z}}_{ciL}{\underset{\_}{I}}_{Bi}{\sinh\left( {{\underset{\_}{\gamma}}_{iL}{d_{LB}\left( {1 - p_{SC}} \right)}l} \right)}}} \right)} \cdot {\mathbb{e}}^{{j\delta}_{B}}}} - {R_{FB}{\sum\limits_{i = 1}^{2}\;\frac{{\underset{\_}{a}}_{Fi}{\underset{\_}{M}}_{i}}{\cosh\left( {{{\underset{\_}{\gamma}}_{iL}\left( {1 - d_{LB}} \right)}\left( {1 - p_{SC}} \right)l} \right)}}} - {{\underset{\_}{a}}_{m\; 0}{d_{LB}\left( {1 - p_{SC}} \right)}{\underset{\_}{Z}}_{0\; m}{{\underset{\_}{I}}_{{Bparal\_}0} \cdot {\mathbb{e}}^{{j\delta}_{B}}}}} = 0$ for double-circuit line where: V _(Bi) denotes the zero, positive and negative sequence of voltage, respectively, measured at the end A, and wherein i=0, 1, 2 denotes the i-th symmetrical component, Z _(ciL) denotes the zero, positive and negative sequence of line impedance, respectively, and wherein i=0, 1, 2 denotes the i-th symmetrical component, I _(Bi) denotes the zero, positive, and negative sequence currents measured at the end A, respectively, and where i=0, 1, 2 denotes the i-th symmetrical component, γ _(iL) denotes the zero, positive, and negative sequence of line admittance, respectively, and wherein i=0, 1, 2 denotes the i-th symmetrical component, R_(FB)—denotes an unknown fault resistance, which a value is received using known iterative method by solving the fault loop formula above, M _(i)—denote numerical coefficients for positive and negative sequences, a ₁, a ₂, a ₀, a _(m0)—denotes weighting coefficients dependent on fault type, a _(F0), a _(F1), a _(F2)—denotes share coefficients dependent on fault type, Z _(0m)—mutual coupling impedance for zero sequence, I_(Bparal) _(—) ₀—denotes zero sequence current from the station B measured in a healthy parallel line, (B3.) then an equivalent impedance of the SC&MOV at fault stage (Z_(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B)) is calculated using post fault measured voltage and currents from the equation: ${\underset{\_}{Z}}_{{SC}\; 1{\_ ph}}^{SUB\_ B} = \frac{\Delta\;{\underset{\_}{V}}_{SC\_ ph}^{SUB\_ B}}{{\underset{\_}{I}}_{SC\_ ph}^{SUB\_ B}}$ and an equivalent impedance of SC&MOV bank at pre-fault (Z_(SC1) _(—) _(pre)) is calculated from the equation: $\mspace{79mu}{Z_{{SC}\; 1{\_ pre}} = \frac{\Delta\;{\underset{\_}{V}}_{{SC}\; 1}^{pre}}{{\underset{\_}{I}}_{{SC}\; 1{\_ B}}^{pre}}}$      where: $\mspace{79mu}{{{\Delta\;{\underset{\_}{V}}_{{SC}\; 1}^{pre}} = {{\underset{\_}{V}}_{{SC}\; 1{\_ B}}^{pre} - {{\underset{\_}{V}}_{{SC}\; 1{\_ A}}^{pre}{\mathbb{e}}^{j\;\delta\; A}}}},\mspace{79mu}{{\underset{\_}{V}}_{{SC}\; 1{\_ A}}^{pre} = {{\cos\;{{h\left( {{\underset{\_}{\gamma}}_{1\; L}p_{SC}l} \right)} \cdot {\underset{\_}{V}}_{A\; 1}^{pre}}} - {{\underset{\_}{Z}}_{c\; 1L}\sin\;{{h\left( {{\underset{\_}{\gamma}}_{1\; L}p_{SC}l} \right)} \cdot {\underset{\_}{I}}_{A\; 1}^{pre}}}}},{{\underset{\_}{V}}_{{SC}\; 1{\_ B}}^{pre} = {{\cos\; h{\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right) \cdot {\underset{\_}{V}}_{B\; 1}^{pre}}} - {{\underset{\_}{Z}}_{c\; 1L}\sin\;{{h\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 1}^{pre}}}}},{{\underset{\_}{I}}_{{SC}\; 1{\_ B}}^{pre} = {{{- \frac{1}{{\underset{\_}{Z}}_{c\; 1\; L}}}\sin\;{{h\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{V}}_{B\; 1}^{pre}}} + {\cos\;{{h\left( {{{\underset{\_}{\gamma}}_{1\; L}\left( {1 - p_{SC}} \right)}l} \right)} \cdot {\underset{\_}{I}}_{B\; 1}^{pre}}}}},}$ Zc1L denotes the characteristic impedance of the line for the positive sequence, γ _(1L) denotes the propagation constant of the line for the positive sequence, in order to determine the final result (d_(B)) which is chosen in the following steps: first, check whether the fault resistances (R_(FB)) is of positive value and if not, then subroutine II is rejected, next, features “real” and “imag” of the estimated equivalent circuit of the SC&MOV is checked for the following dependence: real(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B))>0, imag(Z_(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B))<0 and imag(Z _(SC1) _(—) _(ph) ^(SUB) ^(—) ^(B))>imag(Z _(SC1) _(—) _(pre)), then it is determined that the fault has occurred between the station (B) and the SC&MOV and the result (d_(B)) indicates the distance to fault, and the information about the real fault case is supplied to fault locator FL.
 2. A protection relay equipped with a fault locator comprising means for performing the steps of claim
 1. 